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Abstract. We propose a system of equations that defines Weierstrass-Jacobi's eta- and 
theta-constant series in a differentially closed way. This system is shown to have a di- 
rect relationship to a little-known dynamical system obtained by Jacobi. The classically 
known differential equations by Darboux— Halphen, Chazy, and Ramanujan are the dif- 
ferential consequences or reductions of these systems. The proposed system is shown to 
admit the Lagrangian, Hamiltonian, and Nambu formulations. We explicitly construct 
a pencil of nonlinear Poisson brackets and complete set of involutive conserved quanti- 
ties. As byproducts of the theory, we exemplify conserved quantities for the Ramamani 
dynamical system and quadratic system of Halphen— Brioschi. 
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1. Introduction 

In this work we propose a description of the classical Jacobi's ^-constants and Weierstrass' 
77- function by means of closed and Lagrangian/Hamiltonian ordinary differential equations 
(ODEs). By simple transformations or reductions these equations lead to many well-known 
differential systems. Among these are the Darboux-Halphen system [11, 16], some its 
modifications [25, 3], the Chazy equation [10], and also a Jacobi system of ODEs [18] which 
has not received mention in the modern literature in the context. For both Jacobi's system 
and equations defining the 77-series we work out the Hamiltonian formalism and show 
that they admit a pencil of (compatible) Poisson structures in the sense of Magri [20] and 
formulation as the generalized Nambu mechanics [22] with a certain 4-bracket. 
The three Jacobi's theta-constants are defined by the classical series 

OC GO OO 

Mr) ■= e ^ iT ]T fe e (fc2+fc) ^ 3 (r) := Mr) := £ fe (-l)*e* 2 ™ 

—00 —00 —00 

and the Weierstrass 77-function is defined by the series 

f 1 00 „2fc7rir 

2 I 1 \ ^ e 



24 (1 - e 2fc7rir-J 

Here, the 'time' r is considered to be a complex variable belonging to the upper half-plane 
H + : 8(t) > 0. These series appear in various problems of theoretical physics because of 
their numerous and deep differential properties [3, 10]. Let us mention some of them. 

Three ^-constant series satisfy the following differential identities for logarithmic deriva- 
tives of their ratios: 

dr tf 3 4 4 ' dr i? 4 4 2 ' dr tf 4 4 3 
Yet another and very well-known identity is the sum of logarithmic derivatives: 

^2 h h = 3i 

1?2 1?3 ^4 * V 

(dot stands for the r-derivative) . If we introduce a notation for these derivatives, say 

(i » 

then the quantities (X, Y, Z) satisfy the 3rd order differential system 

X = (Y + Z)X -YZ, Y = (X + Z)Y-XZ, Z = {X + Y)Z-XY, (2) 

which is widely known as the Halphen system [16, p. 330-331]. This system is frequently 
named as the Darboux-Halphcn system though Darboux himself wrote down only differen- 
tials [11, p. 149]: 

C{dA + dB) = B{dA + dC) =A(dB + dC). 
These can be written in the form 

dA _ dB _ dC 

A(B + C) -BC ~ B(A + C) - AC _ C(A + B) - AB _ t 

which is equivalent to the system (2). 

Remarkable applications of Eqs. (2) were initiated in the 1990's by M. Ablowitz ct all 
[7, 1] in connection with reductions of the self-dual Yang-Mills equations. These equations 
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usually provide the main physical motivation for studying both the ij, -^-series and allied 
modular objects. However, applications go beyond the Yang-Mills theory. In succeeding 
years the system appeared in the vacuum cosmological Bianchi-IX model [10, p. 143, 147], 
[3, p. 577], [1], theory of 2-monopolc moduli spaces [5], and many other areas of mathemat- 
ical physics [19]. System (2) has also varieties. One of them is the Weierstrass-Halphen 
dynamical system for Weierstrass' invariants g 2 , g 3 , and ?/-series: 



-p = -(8 fl2 ,-12 fl3 ), -§--{^--giy T ^-{^-- 6 9 2 )- (3) 
It is known that these invariants 

aw : = E (2mr U )4 • = 140 E c^W «»> m ^ * (°- °» 

n.m v 7 n,m v y 

are related to the ^-series by the standard polynomial formulae 



5 3 (r) - ^{^W + ^(r)}{^(r) + <(r)}{^(r) - 0*(t)} 



(4) 



and the series themselves satisfy the well-known Jacobi identity 

^(t)=^(t)+<(t). (5) 

In different notation and (number-theoretic) definition for function series, system (3) is 
known as the Ramanujan system of differential equations [25, 9] for modular forms 

12 12 21fi 

E 2 (r) = - V (r), E 4 (t) = - 92 (t) 7 E 6 {t) = — g 3 (r) . (6) 

TT Z 7T 4 7T° 

Ramanujan's system is sometimes referred to as the Eisenstein system of differential equa- 
tions [9], though Eisenstein himself had not derived it [12]. Further discussions of and 
bibliography to the systems mentioned above can be found in [1, 3, 4, 7, 10, 13] and refer- 
ences therein. 

1.1. Motivation for the work. Dynamical variables for all the systems above are ratio- 
nally expressed through the ^-variables. Therefore, the inverse transformations will involve 
the multi-valued functions, as further examples show. Recently Ablowitz, Chakravarty, 
and Hahn [4] called attention to yet another instance which is more interesting and comes 
from the equations for modular forms on group To (2). This is the Ramamani system [24] 
(Sect. 5.1.1); it was also considered in works [21, 13]. In this case, the relation between dy- 
namical variables and the r\- variables is not obvious because it is given by a duplication 
of the T-argument in forms (6) [4, 13]. If we make use of the duplication rules 

»7(2r) = \v{r) + ^{^W + ^(r)} , g 2 (2r) = -\g 2 {r) + ^{4(r) + ^(r)} 2 (7) 

we arrive again at substitutions of a rational type (see Sect. 5.1.1 for further details). 
Though these rules have not appeared explicitly in the literature known to us, they can be 
established by standard techniques. In a g-series notation these identities can be found in 
Ramanujan's notebooks and (in number-theoretic notation) implicitly have been tabulated 
in [21, Table 1]. Apart from inversions of that kind substitutions one should mention the 
fact that 3-dimcnsional systems, e.g., (3), present generically subsystems or reductions of 
the 4-dimcnsional ones because differentiations intertwine equally all the four objects 
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and 77. In particular, we display here a version of relations which, besides their symmetrical 
form, close differentially these objects. 

Proposition 1. The canonical Jacobi's 'd-constant series satisfy the closed differential iden- 
tities upon adjoining the Weierstrass rj-series: 

(8) 

Derivation of these formulae uses computation of the theta-derivatives through the 
derivatives of g 2 , 53, that is Eqs. (3) and (4), followed by applying the symmetrical identity 
(5). 

To all appearances, these identities, direct consequences of standard relations as they 
are, do not appear explicitly in so extensive literature on thcta-functions though we found 
first three of them in Appendix A to monograph [19]. 

On the other hand, considering (8) as a dynamical system, its integration, as we shall 
see, encounters serious difficulties. Moreover, system (2), under the definition (1), is not a 
consequence of Eqs. (8) but holds only upon restricting to the constant level surface (5). 
To put it differently, the modes of embedding functions 1?, rj into differential systems are 
not unique and Eqs. (8) require modifications. All this will be the subject matter of further 
consideration. 

For the reasons given above it is essential to have a comprehensive description for dif- 
ferential properties of the canonical ?y-series as such. In particular, it is of interest to 
find a Lagrangian and Hamiltonian formulation for these systems. This would provide fresh 
insight into properties of the theta-constants. This circle of questions, as applied to an 
equivalent of the system (3), is addressed in the work [9] by D. Chudnovsky & G. Chud- 
novsky and, to the best of our knowledge, this is the first paper 1 where the problem of 
Hamiltonian treatment for dynamical systems of modular type was raised. These authors 
proposed a 4th order differential system [9] and its reduction to equations of the order 3. 
Although their Hamilton function is rather ingenious and correct, the proposed reduction 
for system (3) is not preserved; in notation of work [9] on p. Ill the reduction is defined by 
the constraint A = 1. In other words, this reduction is satisfied only by a trivial solution. 
It should also be noticed that system (3) and its varieties have a nice interpretation as a 
halfway between integrable systems and the ones with a chaotic behavior. However, we 
do not touch on this kind of problems here because questions of dynamics and transitions 
from 'exactly solvable but not completely integrable' flows to the ergodic ones are the main 
subject of the work [9]. In the same place the extended bibliography is given. 

1.2. The paper content. The subsequent material is organized as follows. In the next sec- 
tion (Sect. 2) we discuss the correlation between identities like (8) and a little-known Jacobi 
dynamical system for Legendre's complete elliptic integrals. We give explanations as to why 
the change from differential //-identities to those that should be thought of as defining 
ODEs is not a trivial question. We write down the simplest version of such ('integrable') 

^See also important comments on p. 5709—5711 of [14] concerning the Darboux-Halphen system (2) and 
its relation to Euler's equations and the Lotka-Voltcrra system; in the same place the detailed comments on 
Poisson structures for the 3-dimcnsional system (2) are presented. In work [8] their construction by means 
of multivalued integrals is discussed and a generalization on Nambu's brackets is proposed. 
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ODEs. Section 3 is technical; it is devoted to explicit integration of system (19) and Ja- 
cobi's system (9). Method of solution invokes standard Legendre's and modular techniques 
and we present results in both of these forms: the 'linear' k- and 'nonlinear' (modular) 
r-representation. In Sect. 4 we explain how these techniques can be exploited in order to 
derive the transcendental multivalued integrals. In Sect. 5 we exhibit explicitly these ob- 
jects for all of the systems mentioned above and use them when constructing Lagrangians 
and the action functional. Complete Hamiltonian formulation to the systems under study 
is expounded in Sect. 6. The found Poisson structures turn out to be non-obvious (none 
are simplectic) and may form compatible pencils; we also describe the genesis of a rational 
degenerate Poisson bracket from a Nambu 4-bracket and possible transitions between var- 
ious Poisson brackets. Section 7 contains a generalization; we complete the theory for the 
Halphen-Brioschi quadratic ODEs. The last section 8 (Appendix) contains some historical 
remarks on Jacobi's system. 



2. ODES DEFINING ^-CONSTANTS 

2.1. On symmetrical system (8). As we mentioned above all the varieties of dynamical 
systems under consideration are algebraically related to each other. In this respect equations 
(8) stand out because this system alone represents the 77, -^-constants. However, point 
transformations between dynamical variables are not unique and resulting ODEs for 77- 
variables may contain parameters. In connection with this ambiguity it is of interest to 
consider an elegant dynamical system which was derived by Jacobi. In Jacobi's record [18] 
it is as follows: 

^ = -166A 2 , 
8h =abA ' 




where h = \ tt\t and the restriction 



0^ = 16(1-26) (10) 



is assumed to be imposed. All the information concerning this system (including solution) 
has been detailed in the next section and Appendix (Sect. 8) contains additional comments 
on original motivation of Jacobi. Jacobi deduced Eqs. (9) as a set of differential identities 
between classical objects of Legendre's 'elliptic theory' [18, 16, 27]: 

1 1 

K(k) = f dX = , K i {k ) = / dX , (ii) 

J v /(l-A 2 )(l-fc 2 A 2 ) J y/(l- A 2 )(A 2 -/fc 2 ) 

k 



m = J \j l -T^dX, E'(k) = J sj 1 f:* )X ' d\. (12) 



A simple computation, based on the i9, 77-representations of objects (11)— (12) appearing in 
Jacobi's definition of variables {A,B,a, b} — this is Eqs. (59) — shows that 

^2 ~i „q4 „q4 „q4 



A = r 3 , b = 



'4w-4 " 1 ^1- &=2 II' (i3) 
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Now, if we drop out the fix constraint (10) and consider (13) just as a point change in 
Eqs. (9), we shall not arrive at symmetrical system (8). We may also insert into the change 
(13) some parameters, say 

and yield different forms to resulting ODEs but we never get the system (8) in this way. 
(Converse is of course also true: symmetrical system (8) does not entail Jacobi's equations 
(9)). Any of such ODEs will be integrable in terms of 77-series since they were obtained 
from (9) by coordinate changes of dynamical (phase) variables. The changes are generally 
algebraic, i. e., multi- valued in both directions. In this respect Jacobi's system (9) is not the 
best variant because choice of the phase variables, in this case, would lead to replacing the 
'simple' K(k) with the 'cumbrous' k(\J '\ ~ \ ^jS+^ ^j ( SCG P ro P os ition 6 further below). 
In other words, the search for a representative defining Jacobi-Weierstrass' series by a 
system of ODE's is not a trivially solvable problem and we need to choose, in some sense, 
'natural and optimal' version for such a system (call it canonical one). It must reflects 
the principal property of the series, namely, the property of being uniformizing for other 
algebraic versions (2), (3), (9), or the like. 

For this purpose, however, symmetrical form (8) is apt to be not a good candidate 
because it is not amenable to integration and we failed to find out its complete integral. 
That such a strong distinction between systems is inherent in the nature of the case (8) 
will be apparent from the consideration of their algebraic integrals as algebraic curves in 
homogeneous coordinates $2 : $3 : $4- 

Proposition 2. The identities (8), being considered as a dynamical system, have an alge- 
braic integral U given by the following rational function of $ 's: 

This integral generalizes Jacobi's identity (5) ifU^=Q. 

Turning now the restriction (10) into algebraic integral (see Eq. (21) in Sect. 3), we 
observe that equation (21), under generalization (14), has genus 9, whereas integral (15) 
is a curve of genus 19. The best we have succeed in solution of system (8) is its partial 
resolution in terms of elliptic functions. In a nutshell, this procedure is as follows. 

Let us change notation U n- U 2 and rewrite integral (15) in form of the elliptic curve 2 

2U 2 xy= (y-x-2) 3 , x = 2^|, V = 2^. 

Hence it follows that the pair (x, y) is parametrized by Weierstrass' (p, p')-functions and 
this curve can be transformed into the canonical Weierstrassian form 

p'(u) 2 =4p 3 (u)- 52 p(u)-<7 3 . 

The computation is rather simple and we obtain 

1 U 2 1 U 2 

x = -p'(u)-p(u) + — -1, y = -p'(u) + p{u)- — + l, (16) 



2 An analogous transformation to the fourth powers of i?'s in integral (21) for Eqs. (9) leads to a zero 
genus curve and no elliptic functions appear in this case (see Sect. 3.2). 
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where constants g 2 , g% are expressed through the integral U : 

J2 12 , W3 21 g g 

Therefore formulae (16) substituted into Eqs. (8) must cause this system to become a t-cvo- 
lution of the uniformizcr u = u(t). Indeed, after some algebra we derive that 

36 U 1 du „ , . tt9 



and therefore 



$ 4 rfr + const . 



12p(u)-£/ 2 36 £/ 
The left hand side of this equation is easily integrated because 

I _^_ = C (u-x ) -C(u + x) + 2C(x), 
where 12p(x) = £/ 2 , p'(x) = ±J7, and C( u ); °"( u ) are the standard Weierstrassian functions 



associated with the basis p(u), p'(u) [16, 27]. We get 
3 J V(u-x) 2CMu 



ln^ -) f e^ wu ^ = mdr + const 

7T1 [cr(u + x) J J 

but integral in the right hand side requires further integration of the system. This last step 
is unknown. 

If equations (8) are indeed non-intcgrable then situation is a manifestation of the mere 
fact that the differential identity for a function and differential equation are not one and 
the same. The function u = — -4- hi(z 2 + -) solves the equation u" = 2u 3 + zu — 2 whose 
general integral is, however, not representable in terms of any known functions or integrals 
of them; this is the 2nd Painleve transcendent [10]. Here is a less trivial example. The 
function u = 4-ln{Ai(z) + aBi(z)} contains, like our ry-solutions, special functions and 
a free constant; functions Ai(z), Bi(z) satisfy the Airy equation ip" = zip. Here, we again 
arrive at the P2-transccndcnt u" = 2u — 2zu + 1. 

2.2. An integrable modification of system (8). Returning to the question of canon- 
ical representative for ODEs defining 77-series, we choose the following modification of 
equations (8): 



dr 



(17) 



Comprehensive explanation as to why the 'defining 1?, 77-equations' should have such a form 
has been detailed in work [6]. System (17) has even an (integrable) extension which is 
described in the same place. It is of interest to observe that all the previous dynamical 
systems contain in effect only squares of ^-constants. For this reason, in the sequel it will 
be convenient to renormalize variables d,ri and adopt the following notation: 
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Then Eqs. (17) acquire the form 

x = (u + y 2 + z 2 )x, z = (u — 2y 2 + z 2 ) z , 

(19) 

<■ , 2 r, 2\ - 24,224 

y = (u + y — 2z )y, u = u— y+yz— z, 

which, along with the Jacobi system (9), will be the main subject of further analysis. Apart 
from simplicity and the symmetry y ±z, there are some additional properties justifying 
the study of canonical system (19). 

First of all, the function u, independently of (x, y, z), satisfies the famous Chazy equation 

ii = 6(2uu — 3 it 2 ) , 

(proof is a direct calculation) which cannot be said of 77-solution to the symmetrical version 
(8). For the latter, the function ^r) solves this equation only if the [/-integral (15) is equal 
to zero (consequence of Proposition 2). Similarly, functions y and z also satisfy a third (not 
fourth) order ODE. This is the known Jacobi C-cquation (60) [18, p. 186]: 

C 4 (ln C 3 C TT ) 2 T = 16C 3 C TT + 36, C=- or -. (20) 

V z 

The above mentioned symmetry involves only functions y and z other than the function 
x. Therefore general solution 1/x(t) docs not satisfy this Jacobi's equation. Equations 
(19) entail that functions x(t) and x{t) times a constant satisfy a common ODE. Hence, 
making the transformation C H> const-C in (20), one infers that the solution 1/x(t) satisfies 
equation (20) wherein 36 should be replaced by a free constant. One easily derives 

C 4 (\nC 3 C TT ) 2 T -l6C 3 C TT = (e^r^V C := - 

\ x 1 J x 

but right hand side of this equation is a constant indeed. Explanation to this fact will be 
apparent from Sect. 3.3 wherein we give a complete integral to the system (19). Thus the 
function C satisfies the 4th order ODE 

(C 4 (lnC* 3 CV r )?) = (l6C* 3 CV r ) 

which is checked by a straightforward substitution. 

3. Explicit solutions and technicalities 

At first, let us integrate Jacobi's system. From (9) it follows that ada = —1696 and this 
equation yields an algebraic integral that replaces Jacobi's restriction (10): 

I 2 = a 2 +326 => 7=0. (21) 

Therefore b is expressed via the function a which in turn satisfies a simple differential 
consequence of (9), namely, the 3rd order equation 



l hhh 



3oL 1 a 2 + 3I 2 



al 2 af 2 (a 2 - I 2 ) 2 ' 
This is a variety of the standard differential equation for Legendre's modulus A := k 2 (r): 

X TTT 3 Arr 1 A 2 — A + 1 

"aT~ 2 Af ~ ~2 A 2 (A-1) 2 ' [ ' 
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It immediately follows that there is bound to be a linear fractional change between variables 
a and A transforming Eq. (22) into equation for a and vice versa. This simple computation 
gives 



21 

Using the well-known ^-constant representation for function X(h) [27], we get 



a = I - 2I ,i 



f(Srf). (23) 

where {a, /?, 7, 5} are free constants with aS ^ /3j. The further integration for the variables 
{A, B} can be continued in two ways. The first one is to make use of rules for differential 
computations (17) of the ?9-series. Applying them to just found expressions for a(h) and 
b(h), we get expressions for A(h), B(h). The second way is to linearize the system because 
any Schwarz's equation is known to be related to a certain linear ODE. We shall give 
solutions both in h- and ^-representations. 

3.1. Associated linear ODEs. Using (23), we have the obvious transformations between 
pairs (a, b) and (k,I): 

a = I-2Ik 2 , 8b = I 2 k 2 (l-k 2 ). (24) 
This allows us to bring (9) into the form 

A = 2A 2 B, B=ll 2 k 2 {l-k 2 )A 3 , k = -Ik{l - k 2 )A 2 , 1 = 0, (25) 
8 2 

where we let the dot above a symbol denote an /i-derivative. We regard this system of 
equations as an intermediate equivalent of Jacobi's system (9) because of its relation to 
linear ODEs. Indeed, as it follows from (25), the quantities A and B, as functions of k, 
satisfy the two linear equations 

dA _ 4 1 n dB - 1 u a 

Ik -Jo^¥jk B > ~dk-\ hA (26) 

and their consequences 

k(k 2 -l)A kk + (3k 2 -l)A k + kA = 0, k{k 2 -l)B kk -(k 2 -l)B k + kB = 0. (27) 

Since k is Legendrc's modulus, it is naturally to expect that these ODEs are intcgrablc in 
terms of functions (11)-(12). 

Proposition 3. Canonical Legendre's elliptic integrals (11)— (12) are differentially closed: 
dK K E dK' kK' E' 



dk ' k (fc 2 -l)fc' dk 1-fc 2 (fc 2 -l)fc' 

dE _ K E dE' kK' kE' 

Ik ~ ~~k + ~k' ~dk ~ 1 - k 2 + k 2 - T 

This system, being considered as a dynamical one, has the general solution 

K = aK(k) - (3K'{k) , A" = -yK(k) + 5K'(k) , 

E = aE{k) + (3 [E'(k) - K'(k)} , E' = SE'(k) + 7 [K(k) - E(k)] 
where {a, (3, 7,(5} are free constants. 



(28) 
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Of course, one should bear in mind that the canonical functions (11)— (12) themselves are 
not independent. Rather they satisfy the Legendre identity 

K{k)E'{k) + K'{k)E(k) - K{k)K'{k) = | Vife, 

which is a particular case of the constant level surfaces for system (28): 

KE' + K'E- KK' = ^ (a5 + (3-y) . 

Curiously, this property and Proposition 3 seems to have not been tabulated in the standard 
texts. The second order differential consequences of this system are known. Both K and 
K' satisfy the same equation 

fc(fc 2 -l)5? + (3^ 2 -l)^ + fc* = => 9={K(k),K'(k)} 
ak^ ak 

and common equation solvable by functions E and E' reads as follows 

k(k 2 -l)^r + (k 2 -l)^--ky = * = {E(k),E'(k)-K'(k)\. 

dk z ak 

The two last linear ODEs are not identical to (27) but search for solutions to Eqs. (26)-(27) 
is not a difficult task. In addition to solution pair (24), we obtain that 

A = Aa.K(k)+A-yK\k), 

(29) 

IB = a[E(k) + (k 2 - l)K(k)] - -y[E'{h) - k 2 K'(k)] 

with some free constants a, 7. Wc can now combine the 'fc-formulae' (29) and /i-time 
dynamics to obtain the complete integral of system (9). 



(30) 



3.2. Solution to the Jacobi system. Let us denote 

ah + (3 
T := -. 

7/1 + 5 

Then, by virtue of (22), 

-fh + 5 K(k) 0§(T) v ' 

Make use of the representation for integrals (11)— (12) through Jacobi's 77, ^-constants. The 
canonical formulae for K and K' are well known: 



'3 V 

One can also show that the second pair {E, E'} has the following modular ^.-representation: 
Modifying these formulae for the general ratio (30), we obtain 

jr(*)-l*Sm - ajrw _ l7jt . W _m_E^SL. (32) 
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Adjust the free integration constants in (29) with those of (30) and (31). Then we may 
write 

r 4~T 



A = ±, 



■{aK(k) -ijK'{k)} 



(33) 



and therefore 

B = ± 



^3 {a[E(k) + (fc 2 - l)K(k)] - ij[E\k) - k 2 K'(k)] } . 



Passing to the 77-representation, we arrive at the final form of the sought-for solution. 
Theorem 4. General solution to the dynamical system of Jacobi (9) has the form 



a = I -21 



0*(T) ' 



b = 



7 2 flj(T)flj(T) 
"8" W?) 



A = ± 



TTi 0§( T ) 



7 -fh + S' 



B = ±1 



7T 3 ( 7 /l + J)l9§(T) 



7T 

12 



^(T)-^(T) + ? /(T) + -17(7/* + <*) 



where {I, a, /3, 7, 5} are /ree constants subjected to normalization a5 — (3j = 1. 

3.3. Solution to system (19). One integral for Eqs. (19) is easily found because x is 
absent in three of these equations. Elimination of u shows that the function 



9 2 

7T 1 = 



(34) 



is a constant on solutions of (19), that is integral. This integral is much simpler than those 
we discussed in Sect. 2.1. As for solutions to system (19), these have the most simple form as 
against the other equations we consider. We shall give these solutions in the next theorem. 
The last fact we should mention here is a point transformation from Jacobi's equations to 
the system (19). The simplest way of getting such a transformation is realized through 
the 'linearizing' systems (25), (26) which can be thought of as intermediate equivalents for 
Jacobi's one (9) or (19). Explanation and details have been given in the previous section. 
From now on we change Jacobi's /i-notation and put 



T 



7T + S 



with normalization aS — /3j = 1. 



Theorem 5. The canonical dynamical system (19) defining $ ^-constants and Jacobi's 
system (9) are equivalent. They are related through the following point transformation 



A 



21 



B = 



1 



■y, 

ll(u + y 2 -2 Z 2), 
V 



12 



2 2 

z V 



2z 2 



b = - 



18 



7T 2 x 4 y 4 



The system (19) has the following general solution: 



7T + S 



6 7T 



2i 77(T) 



(35) 



6 7T + 5 



7r ( 7 r + 5) 2 7r + 5 
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where e ^ is the fourth free constant. With e — 0, the solution decouples into the two 
parametric elementary one: 

±1 1 7 2 t + 7<5-1 



x = 0, y 



Proof. The most convenient way to get the point transformation is to exploit the known 
general solution of Jacobi's C-equation (20) [18, p. 186]: 

c -\ = /Tri ^K^T+f) 
V 6 7 r + 5 

and to pass to the intermediate set of ('linear') variables (A, B,k, I) followed by use of 
identities (24), formulae (31)-(33), and elimination of integration constants {a,/3,-f,5} ap- 
pearing in the general solution given by Theorem 4. Omitting the computation details, we 
derive the transformation {A, B, k, 1} — > {x, y, z, u}: 

x=-^kA, y = (l+i)IA, z 1 = 2iI 2 (l~k 2 )A 2 , u = 2A{B~\I 2 (2k 2 -\) A] , (36) 

where 12iJ 2 = /. This change turns (19) into the system (25). Inverting this change, we 
obtain 



1 — i 1 + i J , , 9 , , 9 1 y 2 — z 2 ,9 z 2 

-^prV, B=— (u + y 2 -2z 2 ), I = — — - — , k 2 - 1 - - 

21 2 y 7r x z y z 



and, subsequently, the substitution (35). Making use of solution given in Theorem 4, we 
get the solution for variables {x, y, z, u}. ■ 

Summarizing, an 'integrable' modification of the change (13) is not obvious a priori; this 
is the change (35). In turn, integral nl 2 x 2 = y 2 — z 2 represents a (corrected) version of 
complicated integral (15) and Jacobi's identity (5) turns into a constant level surface in the 
phase space (x, y, z, u). 

Remark 1. From the preceding, incidentally, it follows that equations of the system (2) 
become now the differential identities for all solutions of Eqs. (17) (proof is a calculation). 
Thus, the Darboux-Halphen system (2) is a su&system for Eqs. (17) but is a reduction for 
system (8). See also the last sentence in Appendix (Sect. 8). 

4. Derivation of integrals 

In order to integrate system (9) we made use of its first integral (21). Having a complete 
solution, we can find the two remaining conserved quantities and the fourth 'integral' cor- 
responds to the time shift /n-» h + e. The simplest way to derive the integrals is to use the 
linear fractional formula (31). Indeed, the /i-derivative of this formula gives the equalities 

d ah + pY 1 (. d K'(ky 



^ h + 5) =[dh^hTs) ={ i dhKM ] = (38) 



and therefore expression 
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must be a perfect square. Upon rooting, we get the h- linear function jh + S with coefficients 
depending on dynamical variables {^4, B, k, 1} {A, B, a, b}. Its /i-derivative 

dh dA 8 V ' dB 2 K ' dk 

yields an /i-indepcndcnt constant 7, that is integral 

••• = Ji(A,B,a,b). 

Doing the same for 



\dhah + [3J \ dhK'(k) i 

we get one more integral J2{A, B,a,b) ~ a. Both of these integrals are independent of 
each other since {a, 7} are independent constants. All the calculus with objects K,k,... 
has been described in the previous section and computations are somewhat lengthy but 
routine. Wc therefore omit them entirely. 

Proposition 6. The Jacobi system (9) has the only algebraic (rational) integral I 2 = 
a 2 + 32b and the two functionally independent transcendental integrals 

Ji = AK(k) ■ B - \E(k) + (k 2 - l)K(k)} ■ AI, 

(39) 

J 2 = 4K'(k)-B+ {E'(k) - k 2 K'(k)}-AI, 
where {I, k}, if required, can be expressed via {a, b} by the inversion of formulae (24): 

I = vV + 326, k 2 



2 2 Va 2 + 326 

Integrals J\, Ji are the multi-valued transcendental functions of dynamical variables {a, 6} 
and the linear ones of {A, B} . 

Curiously, the 'very simple' monomial dynamical system (9) has rather complicated 
transcendently algebraic integrals. Another way of derivation of integrals exploits the linear 
equations (26)-(27) and the well-known Wronskian relation for 2nd order linear ODEs. For 
example, the A-equation in (27) has K(k) as its particular solution. Therefore 

/ 7 \ dA dK(k) A 

K(k) ■ — -f-i ■ A \(k 2 - 1 )k = const. 

dk dk J 

Replacing here -^A by B through (26) and using rules (28), we arrive again at the integral 
J\. The choice of K'(k) for a particular solution produces the second integral J2 in (39). 

5. Integrals and Lagrangian 

5.1. Conserved quantities. Lagrangians, Hamiltonians, and Poisson structures for dy- 
namical systems are known to be closely related to integrals of the corresponding ODEs. 
The Hamiltonian formalism for the systems (9) and (19), which we are about to give is 
based on construction of conserved quantities and we first tabulate the complete set of such 
objects associated with equations (19). 

Proposition 7. The system (19) has the only algebraic (rational) integral (34) and the 
two transcendental multi-valued ones 

Ji = - y i u - 2 y 2 + z2 ) K {^) + 3 ^(5) . j2 = ~( u+ y 2+z2 ) K '{i) - s y E '{i) > ( 4 °) 
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that is J i = Ji = 0. The integrals satisfy the identity 

3 



Proof. Straightforward computation with use of Propositions 3, 6 and system (19) itself. 



2 y 



It should be emphasized here that integrability of any modular system is always associ- 
ated with linear ODEs of Fuchsian class; in particular, with the hypergeometric equations 
[23, 3]. This makes inevitable the appearance of transcendentally multi- valued functions like 
K, K' , E, E' . Transition between such a 'linear' and 'modularly nonlinear' r- representation 
was explained in Sect. 3.1. 

Insomuch as the famous system (2) has an extensive literature and applications, it is 
not without of interest to translate just obtained integrals into conserved quantities for this 
system. The authors of work [14] note properly on p. 5709: 'These conserved quantities 
have not appeared in the literature for over a century even though a great deal of works has 
been done in related areas'. These integrals (including a generalization of (2)) are discussed 
in Ref. [8] and implicitly represented there in a form of 'nonalgebraic, transcendental trans- 
formation . . . between' [8, p. 1755] dynamical variables and integrals through solutions of 
a hypergeometric equation and its derivatives. 

Proposition 8. The following is a complete set of (two) conserved quantities for the 
Darboux-Halphen system (2): 



Z 



' X 


- Y 




- Z 


' X 


- Y 




- Z 



X 



'X 


- Y 




- Z 


' X 


- Y 




- Z 



J 2 = -^=K'[ 

Proof. Taking Remark 1 into account and using the definitions (1) and (18), we may write 
system (19) as follows 

X = u + y 2 + z 2 , Y = u + y 2 -2z 2 , Z = u-2y 2 + z 2 . (41) 

Hence 

u=±(X + Y + Z), y 2 = l -(X-Z), z 2 = l -(X-Y). 

Substituting this into (40), we get the statement of the proposition. ■ 

Not so simple but straightforward computation leads transcendental integrals for dy- 
namical system (3). To do this, one uses the standard notation Pjf , Q£(z) for Legendre's 
functions with indices (i/, /i) = (5,3)- Recall that both of these functions satisfy the equa- 
tion [27] 

(1 - z 2 )^' - 2z^ + [v(v + 1) - Y~^}iJ 
Then one can derive and check the following statement. 
Proposition 9. The following expressions 

Ji = -i/g^w {P>?(g 3 w) - (.93 - j7 1 g 2 )w P£(ff 3 tu)} , 

■h = ~^w-{Q l i(g 3 w) - (53 - i^ss)*"^!/ (&«>)}. 

provide the two independent transcendental integrals for the Weierstrass system (3). 



2 

0. 



_L 3 

.93 27 92 
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It is worth noting here that the Halphen system (2) and Weierstrass' equations (3) are 
equivalent by means of the following transformation (A, Y, Z) <=^ (77,52,53) [16, P- 331]: 

6i?7= tt{X + Y + Z), 
-3.g 2 = tt 2 (A 2 +Y 2 + Z 2 - XY - XZ - YZ) , 
-54i 53 = tt 3 (2A -Y - Z)(2Y - X - Z)(2Z - X - Y) . 

As already noted, all these systems are rationally rcprcsentable through the 77-constants, 
i.e., the phase variables of system (17). Hence, in addition to (1), that is the point change 
(41), we obtain a 'correction' of definition (4); it determines equations (3) as a subsystem 
of (17): 

92 = ^{4~ + *2} > 5 3 = ^{2^ - + *$}{20i - <%} (42) 

(no ^2 here). These changes entail an equivalence of integrals: the {Pjf : Q%, PJt, Q-v}- 
objects can be written in terms of {K, E, K', E'} (and vice versa) by the formulae given 
above. Of course, this property manifests itself in the fact that integrals of such a kind may 
be rewritten solely in terms of hypcrgcometric 2 -Fi-functions; we comment this point more 
fully in Sect. 7. 

5.1.1. On the Ramamani system. By way of illustration to the theory we can also consider 
Ramamani's system mentioned in Introduction. This system is satisfied by certain func- 
tion series V(t), V(t), and <2(r) (see formula (3.3) in [4] and formulae (7)-(9) in [13]). 
Translating these definitions into our notation (77,52,53), we obtain: 

V (t) = ±{ 4 7 7 (2r)-r / (r)}, V(r) = ^{2^)-^)}, Q(r) = jL{l6 5a (2r)-0 2 (r)}. 

Using the duplication rules (7), notation (18), and definitions (42), we may, as pointed out 
above, turn these formulae into the rational point transformation. One obtains 

TriP = 2(u + y 2 + z 2 ), niV = 3 (y 2 + z 2 ) , tt 2 Q = 36y 2 z 2 (43) 

and this substitution brings the main system (19) into the system 

? = kP 2 -S), ^=iri(VV-Q), ^ = 27ri(P-V)Q. (44) 
or 2 cIt dr 

This is the Ramamani dynamical system [24, p. 116], [4, (1.8)], [13, (10)]. Its theory, includ- 
ing the search for conserved quantities Ji,J;i(V t V,Q) (exercise), results from Proposition 8 
and an equivalent of substitution (43) written in terms of Darboux-Halphen variables reads 
as follows: 

niV = 2X, iriT = 2X -Y - Z, ir 2 Q = A(X - Z)(X - Y) . 
See Refs. [4, 21] for further information about system (44). 

5.2. Action and Lagrangians. Let us introduce the collective notation for the phase- 
space coordinates: X := (A, B, a, fe) T for system (9) or X := (x, y, z, it) T for (19). We 
are looking for the action functional 

S= [c(X,X)dT (45) 
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in the first-order formalism. Then the most general non-singular Lagrangian has the form 

£(X,X) = Qk (X)X k -H(X), (46) 

where T-L{X) is a Hamiltonian and Q k {X) define the symplectic potential g = g k (X)dX k . 
As usual, varying (45), we get the Hamiltonian equations 

Xn = nnk{x) _n_ ^ nUx) ± n = _J±_ ^ (9), (19), 

where Q = f2 _1 is the Poisson bi-vector dual to the symplectic 2-form 

dg n (X) d Qk {X) 

kn ~ ax k dx™ ■ 

There is a great deal of ambiguity concerning the choice of Lagrangians for a given 
system of equations. On the other hand, no general method is known for constructing 
Lagrangians starting from equations of motion (the so-called inverse problem of calculus of 
variations). Our situation is, however, somewhat special as the systems under consideration 
are integrable (in the sense that was explained in sects. 3.2 and 3.3). The latter fact allows 
us to write down the following ansatz for the Lagrangian: 

c(x,x) = (Xf-i)n + hi 2 . (47) 

Besides the Hamiltonian 'H(X), it involves two additional independent integrals of motion 
Ij = Ij(X) and the quantity M = N(X) obeying condition N = 1. Clearly, M is a linear 
function of r modulo integrals of motion, that is Af(X) = r + const (%, I\, 12). We have 
constructed such a function in Sect. 4. A computation, based on (39) followed by use of 
(33), (25), and (28), shows that 

AJi dh 
where J\ has been defined in (39). The quantity M{X) thus becomes 

N{X) = ~^JT = (u - 2y* + Z 2)K(\) + WEty 

The Lagrangian density is determined up to a total r-derivative and therefore its choice is 
always accompanied by some heuristic arguments (simplicity of Lagrangians, brackets, etc.) 
When deriving the objects above we made use of a 'linearized' equivalent to systems (9), 
(19), and equation (25). Therefore we present result in terms of 'mixed' phase variables. 
The most compact Lagrange function we have found is given by the following statement. 

Theorem 10. The systems (9), (19), and (25) are the Euler-Lagrange equations for the 
following Lagrangian C: 

d ( B 



-^^- S {^ + ^^}-t + {* + H W -« J0 }./- J f.(«, 

Here, we omitted indication of argument in Legendre's integral K(k) and expressions for 
J\ 1 J2{A 1 B 1 k,I) are taken from Proposition 6. Transitions between variables are described 
by substitutions (35) -(37). 
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6. Poisson structures 

The following statement characterizes a non-trivial property of the dynamical systems under 
consideration. 

Theorem 11. Whatever the Hamilton function H(X) is chosen (single- or multi-valued 
analytic function), none of systems (8), (9), or (17), (19) does admit a constant non- 
degenerate Poisson bracket Q. 

Proof. Denote by X the phase-space coordinate vector for any of these systems: X^ — 
V^(X). Assuming the availability of the form FIX = VH(X) with a constant matrix 
we apply intcgrability condition to equations VH = fiX, considered now as equations for 
the Hamiltonian %: 

w k u = n kj v j => w n (n kj v j ) = w k (n nj v j ) => n kj w n v j = n nj w k v j . (49) 

It follows that Jl ■ d x V must be a symmetric matrix for all X. Straightforward computations 
show that this property is compatible with vector fields V's defining the systems (8), (9), 
and (17), (19) if and only if O = 0. ■ 

This proof gives in fact a criteria for availability of a canonical symplectic form (given 
coordinates) and absence of such a bracket suggests to look for non-canonical one. Insomuch 
as (34) is the only single- valued function integral, we have to take it (or function of it) as a 
Hamiltonian. Furthermore, equations of motion do not depend on choice of the Lagrangian 
C but bracket Q does; even though the Hamilton function T-L(X) and coordinates X have 
been fixed. We thus have to choose in (47), except for H(X), the two independent integrals 
Ii,l2(X) in order that the bracket ft be simplest. We put 

£ = (ri-l)n + (\J 1 )- 1 j 2 , (50) 

where A is an arbitrary constant. Formulae of the previous section contain all what we need 
for computation of the bracket f2. 

Lemma 12. Having integrals of motion H, I k , and the object M, the Poisson bi-vector O 
is calculated by the following computational rule: 

VL = M M T , M kn = W k H ■ V n N + Vfcii • V„I a • 

Now, we insert here 

h := (XJJ- 1 , h-=J2 (51) 

and use Proposition 3. 

Theorem 13. Denote X :— (x, y, z, it) T . Then 
(1) Dynamical system (19) is Hamiltonian: 

1 y 2 - z 2 
X = uVU, H = - - = , 



where the degenerate rational (single-valued) Poisson bracket is as follows 

( (u+y 2 -2z 2 )y (u-2y 2 +z 2 )z u 2 -y 4 +y 2 z 2 -z i \ 

-(u+y 2 -2z 2 )y 

-(u-2y 2 +z 2 )z 

-u 2 +y 4 -y 2 z 2 +z 4 



x 

2H 
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(2) Non-degenerate but transcendental multi-valued extension of the uj is given by the 
bracket 



Q = uj + \uj (detfi^0) : 



where 



2 9 

-A 2 

7T 



x 2 

— z xz x 



Mi \ 



y 

-z 2 -(y 2 ~z 2 ) -M 2 
y y y 

-xz -(z 2 -y 2 ) zM 3 

V 1 

-xMi — M 2 -zM 3 



V y I 

and 

Ah := 3y 2 {ER- 1 - l) 2 - z 2 , M 3 := y 2 (3A 2 A~ 2 - 1) + z 2 , 
M 2 := iy^iER- 1 - l) 2 + y 2 z 2 (6EK~ 1 - 5) + 2z 4 . 
(3) T/ie matrix uj is a bracket as well with detu) = 0. T/ie brackets u), uj are compatible 
to each other and have the following Casimir's functions: 

cjVJi = cjVJ 2 = 0, QVH = ujVj\f=0. 
The system may thus be treated as bi-Hamiltonian in the sense of Magri [20]. 

Incidentally it should be observed that degenerate but well-defined rational bracket uj is 
obtained from non-degenerate but multi- valued bracket f2 by a passage to the limit A — > in 
transcendental part of the f2. This procedure can be interpreted as a formal separability of 
canonically conjugated pairs (H,Af) and (Ji, J 2) in Lagrangian (50). Their commutation 
relations (algebra of integrals) are standard: 

{% J 1 } a = {H, J 2 } Q = 0, {J 2 , (AJi)- 1 },, = 1. 

Remark 2. An explicit analog of Theorem 13 for Jacobi's system (9) is obtained with avail 
of transformation law for tensor Q(x,y,z,u) t- > Q(A, B, a, b) under the coordinate change 
A" := (x,y,z,u) T h-> (A, B,a, 6) T =: Y defined by Theorem 5. Coordinate form of the 
transformations reads 

dY J dY p ~ dY k 



dX n dX m y ' ' dX n 

and implies equations 

Y j = W p (Y)— ^ (9). 

We do not display here the explicit formulae since we were unable to find the compact form 
to them. 

It is interesting to note that in addition to the algebraic integral (21) and the rational 
(but degenerate) Poisson bi-vector, systems (9) and (19) admit a symmetry given by the 
linear vector field 

G = 2xd x = Ad A - Bd B - 2ad a - 4bd b . 
This vector field is of course non-Hamiltonian for otherwise it would be generated by a new 
integral. The absence of rational integrals of motion other than / implies that the result of 
action of G on / should be a function of I. Indeed, one can check that 

GI = -41. 
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Also, once the Hamiltonian form for a dynamical system has been found, we can determine 
its invariant volume form V = VdctJl d 4 X . Calculating determinant of the matrix f2, we 
obtain that 

dctn = -^\ 2 Jt ■ x 6 y 2 z 2 (52) 
and, as in the case of Halphen's system (2) [14, 8], the volume form is a polynomial function: 

dxdydzdu = V = —5- dAdBdadb. 



x 3 yz A 2 b 

Clearly, the invariant volume is not unique as one is free to multiply it on any positive 
function of integrals. 

Let us also comment on the relationship of the degenerate Poisson structure lo to the 
Nambu structure 3 . The general Nambu 4-bracket in a four dimensional space reads [22, 26] 

{/i./a./a./J = S- 1 (X) £ ^™V J / 1 V fc / 2 V,/ 3 V„/ 4 , (53) 
where multiplier B.(X) transforms as a density and e is the Levi-Civita symbol with e 1234 = 
1. Setting S = vdet fl times a function of integrals, one can see that the rational Poisson 
bracket u> appearing in Theorem 13 is a reduction of the Nambu 4-bracket with respect to 
the pair of transcendental integrals I\ , I2 '■ 

{fuh} u = {fi,f2,h,h}- 

Now the dynamical system (19) can be viewed as a generalized Nambu mechanics with the 
4-bracket (53) and the triple of Hamiltonians H, I\, I 2 (two of which are transcendental): 

X = {X,H,h,h} & X = {X,U} U - 

More explicitly, 



VdctOe^™ • W e I 2 ■ V„7i 



and expression (52) leads to a polynomial character of the Nambu bracket: 

where the integrals Ik and Jk have been determined in Eqs. (40) and (51). 

We conclude the section with general remarks concerning other non-constant brackets. 
All of them are obtainable from each other by general transformation of the quantities 
appearing in Lagrangian (50): 

N ^ N+F^H, J lf J 2 ), H ^ F 2 (H, J u J 2 ) , J 2 ^ F 3 (H, J U J 2 ),J 2 ^ F 4 {H, J U J 2 ) (54) 

(re- normalization of integration constants). This defines a function freedom of the three 
variables (a,/3,7) ~ (H, Ji, J 2 ). On the other hand, all the dependencies Q(X), including 
possible change of the Hamilton function Ji, are determined by the following modification 
of the line (49): 

v n (n w v j ) - v fc (n^ v 3 ) (v„ n w - v fe n nj ) v j = n nj w{ -n kJ wz, (55) 

where the tensor field 

dV j 

Wl := — — 

k dx k 



3 The comment has been added following a suggestion of the anonymous referee who we wish to thank 
for that. 
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can be thought of as given. Equations (55) are a set of partial differential equations for 
fi(Jf)'s but, thanks to function freedom mentioned above, we may pass from old set of 
variables, say (x, y, z, u), to the new one (N, a, /?, 7) and thereby turn these equations into 
ordinary differential equations in variable N. 

Theorem 14. Denote fi(N; a, (3, 7) := Q(x, y, z, u) and matrix W = W(N; a, (3, 7): 



dX k 



X=X(N;a,P,i) 

where (a, /3, 7) are seen as parameters. Then all the brackets Q(X) = il(N) satisfy the 
linear matrix dynamical system 

^ = WO+nW T (56) 
supplemented with the arbitrary initial condition (bracket) J?(0) = A(a, /?, 7). 

Proof. With use of antisymmetry Qkj = —tyk and Jacobi's identity V„fifcj + Vfcfij„ + 
Vj fl n k = equations (55) may be rewritten as Vj £l n k = — fl n j Wl . Hence 

tl = — Q,W — W T fl and, subsequently, (56) since a = $ = 7 = and fl f2 = — flf2. ■ 

By this means function freedom (54) with the three functions of three variables a = 
H(X), fi = Ji(X), 7 = J2(X) is converted into the coefhcients of dynamical system (56) 
(matrix W) and its initial condition 17(0). As the latter one may take any particular bracket; 
for example, the bracket from Theorem 13. 

7. A GENERALIZATION 

Outlined receipt of derivation of integrals and the 'linear objects' like Af is directly 
extended to more general (Halphen, Brioschi (1881)) quadratic homogenous systems [15] 

i = x 2 +E, y = y 2 + E, z = z 2 + 'E,, (57) 

S := a(y - x) 2 + b(z - x) 2 + c(z ~ y) 2 , 

associated with a hypergeometric equation of the general type 

s (s - 1) *" + {(a + b + 1) s - c} * + ab* = , (58) 

where prime stands for the s-derivative. Integrability of this system (and its generalizations) 
in terms of associated linear equations was considered and established independently in the 
1990s by many authors: Ablowitz et all [3], Ohyama [23], Harnad & MacKay [17]; see 
also Refs. [8, 2]. Parameters (a, b, c) are computed via the hypergeometric ones (a,b,c) 
(correcting a typo in formula (3.6) of Ref. [23]): 

4a = ac + bc-2ab-c, 4b = a 2 + b 2 - ac - bc + c- 1 , 4c = c 2 + 2ab - ac - be - c 

and base definitions for variables (x, y, z) and relations between them and the quantities t, 
s, and ^ read as follows 

_*(£) . c , _ 1)a+b - c+ i^ 2 

x = 1 A ln ^2 (s _ 1)a+b -c + l s c t _ 1)Q+fc _ c+1 ^ 2 j 

2 dr 2 ds I J 
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1 d_ 
2dr 
ld_ 

2 dr 



ln^ 2 (.s - 1) 



a+b-c+1 g c-2 



r.C-1 



(8-1) 



a+b-c+1^,2 



ln# 2 (s - l)«+6-c-l a c =x _ s c( s _ y 



a+b-cq,2 



See works [23, 17, 8] for details. From the formulae above it follows that 

z — x 



z-y 



(x-y)s, 



2 Fi(a, 6; c|s). 



In the framework of these definitions we obtain system (57) and can deduce its integrals. 
Indeed, passing to the general solution s = we have, instead of (38), 



(7r + sy 



- (x-y)s- 



dr *(s) 



V]/2 



Hence 



( 7 t + <S)~^ K ~7=^= 2F 1 (a,b;c\s). 



Take the T-derivative of this expression and make use of the fact that derivative of a hy- 
pergeometric series is another hypergeometric series [27]: 

^{2*1 («> b- c\s)} =^. 2 F 1 (a + l,b + l;c+ l\s) . 

We thus obtain the first integral J\ ~ 7 for Eqs. (57): 

2 — x 



Ji = C7A- 2 Fi o,6;c 



2 - y 



) +CB- 2 F 1 ( 



CB- 2 Fx o+l,6+l;c+l 



z-y 



where 



A= (a + b — l)x — cy — (a + 6 — c+l)z, 



D 



ab (x — y)(z — x) 
c z-y 

i(c-l) 



C = (y- xY {a+h - c \z - yT i(a+b \z - x) 

Assume now that the second (linearly independent of solution to (58) has no logarithmic 
behavior in the vicinity of point s = 0; otherwise we can reorder variables (a;, y, z) with the 
help of the linear transformation s i— > 1 — s or s 1— > s _1 . If the logarithm presents at each 
of the points s = {0, 1, 00}, we fall into Proposition 8. This case corresponds to parameters 
(a, b, c) = 5, l) an d is equivalent to system (2) up to a simple linear transformation [23] 
of the triples (x,y,z) (X,Y,Z). Then we may take the following form for the second 
solution to (58) [27]: 

§ = s 1 - ^ - l) c - a - b ■ 2 Fi(l - a, 1 - b; 2 - c\s) . 
By repeating the arguments above we obtain the second integral for Eqs. (57): 



J 2 = CA ■ 2 Fi 1 - a, 1 - b; 2 - c 



z-y 



where 



A=A + 2(z + y), 



D 



CB ■ 2 F 1 (2-a,2-b;3-c 



C 



,(o-l)(6-l) (x-y)(z-x) 



z-y 



c 



c-2 z-y (z-y) 2 

This completes an integration procedure considered in Refs. [2, 3, 8, 17, 23]. 
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8. Appendix: The Jacobi system 

Since the late 1850's C. Borchardt, being the Editor-in-Chief of Crelle's Journal, began to 
edit and publish the manuscript material kept after Jacobi's death in 1851. In particular, 
in 1857 he published calculations [18, p. 383-394] where Jacobi constructed power series 
developments for his 0(z|r)-functions. The power 0-series are of interest in their own rights 
but not a less remarkable fact is that they produce the nice dynamical systems intcgrable 
in terms of ^-constants. 

Jacobi introduces the four variables (we keep completely to Jacobi's notation in [18, 
p. 386]) 

A=— , B = — -k' 2 —, a = 4(l-2fc 2 ), b = 2k 2 k' 2 , (59) 

TT TT TT 

and shows that these satisfy the dynamical system (9); in doing so Jacobi imposes the 
condition (10) which is of course an equivalent of the relation k 2 + k' 2 = 1 or, which is the 
same, the z?- identity (5). Halphen does not mention system (9) and, to all appearances, it 
has not received mention in the later literature on theta- functions. Jacobi does not restrict 
his consideration to variables (59) and exhibits what is called presently the canonical trans- 
formations, i.e., transformations of dynamical variables preserving the form of equations. 
Here are his versions of the transformations [18, p. 387]: 

2kK 1 2E 4(1 + k' 2 ) u _ 2k' 2 

A — , B — — • , 



k tt k 2 k 4 

2k' K 1 (2E 2K\ 4(1 + k 2 ) , 2k 2 
, B = - , o= V ... ' b = 



k' \ TT TT J ' k' 2 ' k' A ' 

Complete set of differential relations between these and auxiliary variables {k, k', K, E} was 
written down by Jacobi earlier [18, p. 176-177]. As in the previous differential systems (2) 
and (3), dynamical variables {A, £?, a, b} are expressed through the 77, ^-series rationally; 
see formulae (13). 

We also note that system (9) is notable for its homogenous monomial structure. Jacobi 
exploits intensively this fact when deriving the power 0-series; the pages 388-391 of his 
Werke [18] contain a lot of useful formulae along these lines. System (9) is not the only 
dynamical system that was derived by Jacobi in connection with ^-functions; see also [18, 
p. 173-190]. Jacobi did not pose a question about integration of (9) as ODEs, however 
earlier, in 1847, he obtained a complete integral for the 3rd order differential equation 

C* 4 (In C 3 C Tr ) 2 = WC 3 C TT - tt 2 (60) 

satisfied by each of the ^-constants: C — t9(t)~ 2 (Jacobi's notation [18, p. 179]). On the 
other hand, this equation must be a certain consequence of equations (8) whose solutions 
are not only the 77-series. Invoking integral (15), we conclude that Jacobi's equation (60) 
is indeed the consequence of equations (8) with the proviso that U = 0. It is also clear 
that this condition is a necessary one in order that the Darboux-Halphen system (2) be a 
consequence of symmetrical identities (8) as well. 
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